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Abstract 

In this paper It Is presented a detailed numerical Investigation of acoustic 
emission signals obtained from test samples of fibreglass reinforced polymeric 
matrix composites, when subjected to tensile and flexural tests. Various fractal 
Indices, characteristic of the signals emitted at the different structural failures of 
the test samples and which satisfy non-stationary distributions, have been deter- 
mined. From the results obtained for these Indices, related to the Hurst analysis, 
detrended fluctuation analysis, minimal cover analysis and to the boxcountlng 
dimension analysis. It has been shown they can discriminate the different failure 
mechanisms and, threfore, they constitute their signature. 
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1. INTRODUCTION 

In a recent paper Ferreira et al [1] discuss the characterization of failure 
mechanisms that occur in fibreglass reinforced polymeric matrix composites 
when subjected to tensile and flexural loads. The characterization was based on 
the analysis of acoustic emission signals emitted by the composite during the 
process of failure, which constitutes one the most important non-destructive 
testing for the detection of structural flaws in composite materials [2-5]. 

The samples studied were manufactured with E-glass fibre roving reinforced 
DER 331 epoxy resin and its preparation and experimental conditions are de- 
scribed in detail in ref. [1]. Besides tensile tests, flexural tests at three- and 
four-points were also applied and four failure modes have been observed, namely, 
matrix cracking, fibre braking, fibre/matrix debonding and delamination. 

The main purpose of the study was to find the signature of these failure 
modes in the acustic emission signals. In order to identify these signatures, the 
signals were studied by using Fourier spectral analysis and wavelet analysis. Al- 
though relevant information has been obtained from these analyses, the authors 
in ref. [1] have not been able to characterize in a clear way the various failure 
mechanisms. 

Therefore, in this paper we readdress the problem by looking at some fractal 
properties of the acoustic emission signals. In particular, we obtain the fractal 
indices related to the Hurst analysis [6], detrended fluctuation analysis [7], min- 
imal cover analysis[8] and to the boxcounting dimension analysis [9], which will 
be used to characterize the different failure modes. 

These types of analysis have been widely used in the study of random non- 
stationary series ranging from seismic [10] and chmate data, [11] to wind speed 
[12] and financial data [13], and in the study of different music genres [14]. Their 
use in the characterization of acoustic signal has been introduced by Duta and 
Barat [15] in the analysis of ultrasonics backscattered signals obtained in the 
study of single crystal and polycrystalline materials. More recently, Matos et 
al. [16] have used this approach to characterize the ultrasonics backscattered 
signals obtained in the study of the cast iron with lamellar, vermicular and 
spheroidal microstructures. 

The study presented in this paper extends the above mentioned analyses to 
a new tj^pe of acoustic signals, namely, the ones obtained in the acoustic emis- 
sion nondestructive testing. The main objective of the work is to show that the 
parameters determined from these analyses can characterize the failure mecha- 
nisms in composite studied. To this aim and in order to establish the parameters 
to be calculated, we present in section 2 a brief review of the numerical analysis 
used in the treatment of the data, and in section 3 we present and discuss the 
results obtained. 
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2. NUMERICAL ANALYSIS 



The numerical treatment of the signals will be performed on data from A- 
scan, which contains the amplitude of the acoustic emission signals as a function 
of time. The parameters to be determined, as pointed out in the introduction, 
will be obtained from the Hurst analysis (R/S analysis) [6], detrended fluctu- 
ation analysis (DFA analysis) [7], minimal cover analysis [8] and boxcounting 
analysis [9]. 

In order to make the paper self-contained and to introduce the notation, 
we will present a brief reA'iew of the these numerical techniques which will be 
used in the analysis of the temporal series. They will be identified as the set of 
random values {yi}, where the label i corresponds to the time variable, which 
satisfy nonstationary distributions. 
2.1 Hurst analysis 

The R/S analysis will provide information on the temporal correlations, on 
various time-scales, of the data. Given the temporal scries {Ui}, with N terms 
(1 ^ z ^ N), we define the average in the interval n as 

1 " 

<y>n= - yZyi' 



and the accumulated deviation from the mean as 



Yij,n) = J2iyi- <y>n), (2) 

where n varies from 2 to A''. 

From these results, we can also define in the interval n the range R{n) of 
the accumulated deviation in the form 



R{n) = max Y{j,n) — min Y{j,n), (3) 



and the standard deviation S{n) as 



Sin) = ^ ^^^y^-<y>-\ 

Finally, we can obtain the rescaled range i?(n)/S'(n)which should satisfy the 
scahng relation 
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where H is the Hurst exponent [6]. 

In the scaUng regime, the previous expression can be written as 



which defines the amplitude Ah- Although this parameter has no universal 
characteristic, as the amplitudes in the scaling laws in critical phenomena [17] 
where they are related to the interactions, they can be used however as an 
additional parameter to characterize the temporal series. 
2.2 Detrended fluctuation analysis 
The DFA analysis [7] aims to study the temporal correlations by eliminating 
the spurious trends in the data which can conduct to misleading results. The 
method consists initially in obtaining a new integrated temporal series {zi}, 
from the original one {yi}, given by 



i=l 

where the average < y > is defined as 



1 ^ 

In the following step the series is divided in time intervals of width n, and 
an order-Z polynomial is fitted in each interval, and we identify the analysis as 
DFA-Z. Then, the detrended variation function of order I in the interval j, 
is obtained by subtracting the local trend contained in the fitted polynomial, 
and is given by 

jn 

^'M)= E (9) 

i=0-l)n+l 

where z\ is the value from the fitted polynomial. 

Finally, we calculate the mean root square fluctuation F\n) 



F\n) 



int[N/n\ 

^ E (10) 



N 



which should scale as 



4 



F'(n) ~n«, 



(11) 



where a is the scaling exponent. 

The dotrondod fluctuation analysis that wc will present will be restricted to 
the linear case, namely, DFA-1. As in the case of the R/S analysis, eq.(ll) can 
be written in the scaling regime as 

F'(n)=A„n", (12) 



which also defines a new characteristic parameter Aa. 
2.3 Minimal cover analysis 

This method has been recently introduced [8], and it relates the minimal 
area necessary to cover a given plane curve, in a specified scale, to a power law 
behaviour. The scale is introduced by dividing the domain of definiton of the 
function in n intervals of width S. In each interval j (1 ^ j ^ n) we can associate 
a rectangle of base 6 and height A{j) defined as 

Aj = max{yi,ie[j,j + S]}-mm{yi,ie[j,j + 6]}, (13) 



such that the minimal area will be given by 



S{6)=J2AjS. (14) 
In the scaling region, S{S) should behave as 

S{6)~6^-^'^, (15) 

where is the minimal cover dimension, which is equal to 1 when the curve 
presents no fractality. We can also define a new exponent fj, given by 

M = - 1> (16) 
which measures the fractality of the curve and satisfies the scaling relation 

V{6) ~ 5-^, (17) 
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where V{d) is the summmation of the heights of the rectangles 



n 



(18) 



The amplitude A^, as in the previous cases, is defined in the expression 



and it also constitutes a new characteristic parameter. 

2.4 Boxcounting analysis 
The boxcounting dimension, which is one of the best known fractal dimen- 
sion [9], is easily defined and obtained numerically. It can be introduced in a 
general d-dimensional euclidean space, where a hyper-volume is embedded, by 
considering the number of hypercubes of side length 6, J\f{S), necessary to cover 
the entire volume. As (5 — »■ 0, A/'((5) satisfies the scaling relation 



where Db is the boxcounting fractal dimension. 

For non-fractal objects, this dimension corresponds to the topological di- 
mension and, in particular, for continuous planar curves .Ds is equal to 1. 
The amplitude As of the scaling relation is, in this case, given by 



V{S) = A^6-^, 



(19) 



M{6) 



(20) 



and it also constitutes a new characteristic parameter. 



(21) 
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3. RESULTS AND DISCUSSIONS 

A detailed description of the experimental setup and samples used for the 
acquisition of the acoustic emission signal from the different tests is presented 
by Ferreira et al. [1]. Besides the tensile test, the samples were submitted 
to three- and four-point fiexural tests. They have identified four basic failure 
modes, namely, matrix cracking, fibre breaking, fibre/matrix debonding and 
dclamination. In the different tests, the failure mechanisms were a result of a 
combination of these failure modes, and they are presented in Table 1. We also 
present in this Table the acronyms for the different specimens. 

In order to reduce the noise in the data, the signals have been processed 
with an adjacent low-pass filter with five points. For each type of specimen the 
tests were carried out in 03 samples, which correspond to the number of signals 
available for each kind of mechanical failure. 

In Figs. 1-4 we present the various analyses made in a given signal from the 
TEM specimen.. These analyses are representative of the results obtained in the 
study of the other signals. In the Hurst analysis. Fig. 1, the crossover from short- 
to long- time correlations is always present. As can be seen in Figs. 2-4, this 
crossover also exists in the box counting analysis and in minimal cover analysis, 
but not in the DFA analysis. It should be noted that the presence existence of 
this crossover on the fractal analysis is characteristic of a multifractal behaviour. 

Besides the eight parameters H, Ah, a, A^, jJ,, A^^, Db, Ab, we can yet define 
an additional one which corresponds to the standard deviation a of the signal. 
This parameter has been recently introduced in the context of the characteri- 
zation of climate of different regions in the United States from the analysis of 
maximum daily temperature time series [11]. 

From what we have presented, this multi-dimensional space parameter can 
be used to discriminate the various types of mechanical failure. As it will be 
shown in the figures relating the characteristic parameters, there is not a unique 
signature, since different combinations of the indices can lead to the identifica- 
tion of the signals. We have restricted our study to subspaces of the parameter 
space, as they can provide the desired signature of the signals. Explicitly we will 
consider the exponents H, a, fi and Db a.s functions of the standard deviation of 
the signal, a, and also as functions of the logarithm of its respective amplitude 
A. 

These functions, which correspond to projections of the points of the pa- 
rameter space in different planes, are shown in Figs. 5-17. Even considering 
that the data for each type of specimen consisted of three samples only, which 
is a poor statistical sampling, we have calculated the standard deviation of the 
variables presented as the error bar on these figures. From the analysis of these 
figures we can verify that the first discrimination attained is the separation of 
the failures caused by traction from the ones caused by fiexion. This discrimina- 
tion is clearly seen in the diagrams Hi x a, H2 x a, Db2 x <^ and ax a which are 
presented in the Figs. 5, 6, 8 and 11, respectively. This separation is the easiest 
to be obtained since the stress distributions in the samples are very different in 
the two cases, and this has strong effect on the acoustic emission signals. 

By starting from any of these diagrams we can obtain the complete discrim- 
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ination of all the failures mechanisms and even distinguish the results obtained 
from the three- and four-point flexural tests. This can be achieved by using 
the complete set of figures, namely, Figs. 5-17. In terms of these diagrams, all 
possible solutions of the problem are presented in the tree type graph shown in 
Fig. 18. In the branches of the tree, we designate the various failure modes and, 
in the nodes, we show the different two-dimensional diagrams which can lead 
to the desired discrimination. As can be seen on the solution tree, we can iden- 
tify unmistakably all the failure modes from the acoustic emission signals and, 
moreover, show that there are multiple paths which lead to the identification 
we are looking for. 
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Table 1. Failure modes and respective mechanical tests for the 
different specimen types. 



Test 


Specimen type 


Failure modes 


Tensile 


Epoxy (TME*) 


Matrix cracking 


Tensile 


Fibre /Epoxy (TLEV*) 


Fibrebreaking 
Matrixcracking 
Fibre / m,atrixdebonding 


Tensile 


Fibre/Epoxy(TTEV*) 


Matrix cracking 


4-point flexural 


Fibre/Epoxy(F41*) 


Fibrebreaking 

Matrixcracking 
Fibre / matrixdebonding 


4-point flexural 


Fibre/Epoxy(F41S*) 


Fibrebreaking 

Matrixcracking 
Fibre / matrixdebonding 


4-point flexural 


Fibre/Epoxy(F42*) 


Fibrebreaking 
Matrixcracking 
Fibre / matrixdebonding 
Delamination 


3-point flexural 


Fibre/Epoxy(F31*) 


Fibrebreaking 
Matrixcracking 
Fibre /matrixdebonding 


3-point flexural 


Fibre/Epoxy(F31S*) 


Fibrebreaking 
Matrixcracking 
Fibre /matrixdebonding 


3-point flexural 


Fibre/Epoxy(F41*) 


Fibrebreaking 
Matrixcracking 
Fibre /matrixdebonding 
Delamination 


*Specimen acronyms. 


F41S and F31S identify the samples with surface treatment of fibres 
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jure captions 

. 1- Hurst analysis of a signal obtained from a TME specimen. 
Hi and H2 are the Hurst exponents associated with 
short- and long-time correlations, respectively. 
. 2- DFA analysis for the signal used in Fig. 1. 
. 3- Minimal cover analysis for the signal used in Fig. 1. 
/ii and /i2 arc the variation indices associated with 
large and small fractal scales, respectively. 
. 4- Box counting analysis for the signal used in Fig. 1. 

Dbi and Db2 are the box counting dimensions associated 
with large and small fractal scales, respectively. 
. 5- Hurst exponent Hi (short-time correlations) as a function 

of the standard deviation a of the signal. 
. 6- Hurst exponent H2 (long-time correlations) as a function 

of the standard deviation a of the signal. 
. 7- Box counting dimension Dbi (large fractal scale) as 

a function of the standard deviation a of the signal. 
. 8- Box counting dimension Db2 (small fractal scale) as 

a function of the standard deviation a of the signal. 
. 9- Variation index /xi (large fractal scale) as a function 

of the standard deviation a of the signal. 
. 10- Variation index /7,2 ( small fractal scale) as a function 

of the standard deviation a of the signal. 
. 11- DFA exponent a as a function of the standard deviation a 
of the signal. 

. 12- Hurst exponent iJi (short-time correlations) as a 

function of the logarithm of the amplitude Ahi ■ 
. 13- Hurst exponent H2 (long-time correlations) as a 

function of the logarithm of the amplitude Ah2 ■ 
. 14- Box counting dimension Dbi (large fractal scale) 

as a function of the logarithm of the amplitude Ajj^^^. 
. 15- Box counting dimension Db2 (small fractal scale) 

as a function of the logarithm of the amplitude Ad^^. 
. 16- Variation index fii (large fractal scale) as a function 

of the logarithm of the amplitude A^i . 
. 17- Variation index /i2 (small fractal scale) as a function 

of the logarithm of the amplitude A^2- 
. 18 - Tree summarizing all possible solutions for the discrimination 

of the various failure modes, which correspond to different paths 
on the complete graph. The modes are shown on the branches, and 
the diagrams discriminating the modes are attached to the nodes. 
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